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ABSTRACT 

We present a morphological analysis of distant field galaxies using the deep ACS images from the public 
parallel NICMOS observations of the Hubble Ultra Deep Field obtained in the F435W (B435), F606W (V 606 ), 
F775W («775) and F850LP (zgso) filters. We morphologically segregate galaxies using a combination of visual 
classification and objective machine based selection. We use the Asymmetry (A) and Central Concentration 
(C) parameters to characterize galaxies up to Z850,ab< 25 mag. We take advantage of the multicolor dataset 
and estimate redshifts for our sample using the Bayesian photometric redshift (BPZ) which enables us to 
investigate the evolution of their morphological demographics with redshift. Using a template fitting model 
and a maximum likelihood approach, we compute the star-formation rate (SFR) for galaxies up to z — 1.3 and 
its contributions from different morphological types. We report that spirals are the main providers to the total 
SFR. The E/SOs contribution flattens out at z — 1 while the Irr/Pec populations continuously rise to match the 
spirals contribution at z — 1.0. We use the /775-Z850 and V606 _ '775 color-magnitude diagrams to constrain the 
galaxies' formation histories and find that E/SOs show both a population of luminous red galaxies in place at 
z ~ 1 .2 and a bluer and fainter population resembling those of Irr/Pec at similar redshifts. 

Subject headings: galaxies: evolution — galaxies:structure — galaxies:formation — galaxies: elliptical and 
lenticular, cD 



1. INTRODUCTION 

The archeological nature of galaxy evolution studies pre- 
vents us from following individual galaxies over time. In con- 
sequence we are only left with snapshots at different lookback 
times from where we attempt, as excavators, to piece together 
their evolutionary histories from the motion of different ob- 
servable quantities and their inter-relationships through dif- 
ferent cosmic times. As morphology correlates with a range 
of physical properties in galaxies, such as mass, luminosities 
and particularly color, this suggest that their appearance must 
embody some important clues about their formation histories. 
Moreover, the growing acceptance of the notion that the mor- 
phological appearance of galaxies may not be a static property 
set at an early stage of formation, makes it crucial to under- 
stand the flow of the morphological mix of galaxies and their 
observables as a function of redshifts. 

The early morphological studies of HST galaxies 
llCowie et alll995llvan den Bergh et alJl996tlVolonteri et al.l 
2000; Reshetni kov et all2003|) set the stagefo r the Advanced 
Camera for Surveys fACS: iFord et all 120021) to effortlessly 
resolve galaxies at ever fainter limits thus enabling the 
study of the morphological demographics of distant galax- 
ies with larger and deeper datasets (see lElmegreen et al.l 
2004a b, 2005, for a list of recent studies). Since its in- 
stallation, ACS has provided a continuous flow of high 
resolution optical imaging of distant field galaxies over 
significantl y wider areas . Large p rograms such as GOODS 
(Giavalisco etal]l2004bl). GEMS JRix et all l200l a nd the 
ACS/GTO dPostman et ail 120051 iFordetalJ 120021) have 
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spawned several studies that address the color evolution of 
galaxi es, the num ber evolution of red objects with redshift 
tBell et al] 2004), t he ma s s ass embly rate for different 
morphologies (iBundv et a D]2Q i Q i 5j) and constrain the ages 
and m asses of early types llTreu et al.ll2005t Evlenantea u et alJ 
2004). 

In this paper we take advantage of public deep multicolor 
ACS parallel observations to study the morphological color 
and star formation rate evolution of galaxies at depth s com- 
parab le to those of the Hubble Deep Fields (Willia m's et alJ 
1996), but probing areas many times larger. Throughout this 
paper we use Hq = 70 km s _1 Mpc _1 and a flat (fij = 0) cos- 
mology with Q,m = 0.3 and Vt\ = 0.7. Magnitudes are given 
in the AB system across the paper. 

2. OBSERVATIONS 

2.1. ACS Multicolor Imaging 

Our analysis is based on the deep ACS parallel fields of the 
Near Infrared Camera and Multi-Object Spectrometer (NIC- 
MOS) imaging of the Hubble Ultra Deep Field (UDF) (GO 
9803, PI: Thompson). The ACS datasets consist of two sep- 
arate fields around the center of the UDF arranged in a mo- 
saiced pattern determined by the NICMOS ob servati ons of 
30" spacing between pointings (Blakeslee et al. 2004). The 
resulting coverage of the parallel ACS fields as well as the 
UDF and NICMOS observations positions are shown in Fig^ 
Each field comprises 9, 9, 18, and 27 orbits in the F435W 
(B435), F606W (Veoe), F775W (/77s) and F850LP (z 850 ) fil- 
ters respectively. The images are publicly available from 
the STScI archi ve and were processed using the ACS/GTO 
pipeline (Apsis; iBlakeslee et afl l2"003a) which includes ob- 
ject detection, photomet ry and cataloging using SExtractor 
(Bertin & Arnouts 1996). Photometry was calibrated in AB 
magnitudes using the zeropoints from Sirianni et al. (2005). 
Photometric redshift s were al s o com puted using the Bayesian 
estimator BPZ from Benitez (2000). Within the deepest re- 
gions of the fields the 10er limited magnitudes were 28.8, 
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FIG. 1 . — The diagram of the relative positions and sizes of the final mo- 
saiced patterns of the ACS parallel observations of the UDF NICMOS and 
the ACS/UDF observations. 



29.0, 28.5 and 27.8 in B435, V606> '775 an d ^850 respectively 
llBouwens et al.ll2004l) in a scale of 0"05 pixel -1 . The final 
ACS area covered was 55.7 arcmin 2 from both fields. For a 
detailed descriptio n of th e dataset and its processing we refer 
to lBlakeslee et al.N2004l) . 

2.2. LDDS2 Magellan Spectra 

Follow up multi-object spectroscopy on one of the parallel 
fields was performed using the low-dispersion survey spec- 
trograph 2 (LDSS2) on Las Campanas Observatory Magellan 
Clay 6.5m Telescope. A complete description of the obser- 
vations and data reduction analys is of the spectroscopic sam- 
ple is given in lMotta et alJ (120051) : here we briefly summarize 
these observations. For Field #1 three masks were observed 
in November 25 and 26, 2003 from which spectra for a to- 
tal of 56 galaxies were secured. From these we secured red- 
shifts for 51, and for a subset of 37 galaxies line strengths and 
equivalent widths (EW) were measured. Targets for the masks 
were uniformly selected from galaxies within the field with 
'775< 23. The spectra were reduced using a combination of a 
customized IDL reduction software techniques optimized for 
the extraction of background limited data ( Frve et al. 2002) 
and IRAF tasks for the wavelength and flux calibrations. The 
final integrated exposure times per mask were between 7200s 
and 4800s at a central wavelength of 5500A and spectral cov- 
erage of 4000-9000A at 5.3A pixel" 1 dispersion. 

2.3. Morphological Classification 

Rather than using the color of galaxies as a surrogate for 
morphological type, in this paper we exploit the superb abil- 
ity of ACS to resolve and morphologically segregate distant 
galaxies with no a priori built-in biases in the selection cri- 
teria. We classified galaxies according to their morphologies 
down to Z85o.ab < 25 using a two-pronged approach based on 
visual inspections and automated classifications. First, galax- 
ies were visually inspected and classified using the Medium 



Deep S urvey (MDS) sche me as described in lAbraham et all 
( 1996 b| and mo dified by Menant eau et alJ i 19991) ("see also 
Treu et al. 2005). Within this system we segregated galaxies 
into three broad categories: i) E/S0 which include E, E/S0 and 
SO; ii) Spirals: including Sa+b, S and Sc+d and z'z'z) Irr/Pec 
including Irregular and Peculiar types. We also limited our 
sample to objects with > 90% confidence in their BPZ red- 
shift estimates. We elaborate on this in the next section. A 
final catalog containing 1228 galaxies was visually inspected 
by one of us (FM) in both ACS fields. 

Secondly, we independently investigated the morpholo- 
gies of galaxies using automated classifications. For this 
we calculated the central concentration (C) and asymmetry 
(A) parameters for a l l galax ies using the definitions from 
Abrah am et alJ ( 11994 fl996bl) . Our analysis is based on the 
well-known relationship between the A and C structural pa- 
rameters and traditional visual morphologies (| Abraham et al. 
1996aB iBrinchmann eUd1fl998l: IConsehcdl2003T) ' We com- 
pute asymmetry and concentration using PyCA, our own au- 
tomated Python software. The code is designed to compute 
A and C from the SExtractor products generated by Apsis, the 
detection catalog and associated segmentation image. How- 
ever, it can be used with any standard SExtractor products and 
it is publicly available 5 . 

Concentration is computed as the ratio between the flux at 
30% of a given radius Q.3R to the total flux up to that radius 
R. Formally it is defined as, 

C = f(03R)/f(R), (1) 

were f(R) = 2tt Jq I(r)dr is the integrated flux within the ra- 
dius R. The determination of the radius R up to where to in- 
tegrate the light is of importance particularly for the asym- 
metry values. As we sample galaxies over ever larger cosmic 
ages, the surface brightness of galaxies at higher redshifts be- 
come heavily affecte d by the well-known (1 +z) 4 cosmologi- 
cal dimming ( Tolma nl 1 934t ISandage & Lubir «200ll) . There- 
fore to consistently measure structural parameters of galax- 
ies over similar physical areas at different redshifts, we chose 
an aperture ra dius defined as R n = 1.5 x r p , with 77 = 0.2, 
where r p is the Petrosian ( 1976) radius and rj(r) = I(r) / (I(r)} 
is the ratio of the galaxy average surface brightness, I(r), 
in an annulus of radius r and the mean value up to the 
same radius, (I(r)) (see iPapovich et alJl2003l iBlanton et alJ 
2001; Menanteau et al. 2004, and references therein). The 
advantage of using a Petrosian radius over traditional sur- 
face brightness limits ones, such as the Holmberg radius (i.e. 
Abrah am et alJl!996at iMenanteau et aljfl999l) is that it only 
depends on the galaxy light profile and i t is therefore in depen- 
dent of the redshift of observation (see Bers hadv et alJl2000l 
for a discussion). Comparatively accurate results can be ob- 
tained using the SExtractor Kron radius, which in most cases 
coincides with R p . 

Similarly, asymmetry is computed as half of the ratio be- 
tween the absolute value of subtraction of the image, 7y, and 
its 180° rotated image around its center, ifp to the original im- 
age, both contained in an elliptical mask of semi-major length 
Rp. To avoid variations arising from noise features, the pivot 
point of the galaxies was computed after the images were 
smoothed with a Gaussian kernel of 1 pixel and then self- 
subtracted. Formally, our working definition for asymmetry 

5 It can be downloaded from http://acs.pha.jhu.edu/~felipe/PyCA 
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FIG. 2. — The asymmetries and central concentrations values keyed to their 
visual morphological classes for all 768 objects in the Field #1. E/SO systems 
are represented by ellipses; Spiral galaxies by spirals and Irr/Pec as asterisks. 
The dashed lines represent the selection limits based on A — C that maximizes 
the recovery of visual classes and minimizes the contamination from other 
morphological types. Numbers inside the figure represent the recovery and 
contamination percentages of visual classes and from other classes respec- 
tively using the A and C based selection. 



over a galaxy with an intensity matrix /;,- can be expressed as, 

A .lEWHSI-». (2) 

Because we use the absolute values of the self-subtracted 
galaxy residuals, small variations in the image background 
can appear even in the most symmetric objects, consequently 
we add a factor b to correct for this artificial A signal intro- 
duced, such as b = v / 2er s kyA'pix, where Af p ; x is the area inside 
the ellipse mask and a^ y is the sky mean standard deviation. 

We compute A and C for all objects in our galaxy sample 
to a limiting magnitude Z85o,ab< 25. To avoid morphologi- 
cal variations resulting from redshifting bandpass (i.e. mor- 
phological Tif-correction), we compute A and C using the 
image in the closest available band to the rest-frame B de- 
rived from the photometric redshift information. In Fig. [2] 
we show the distribution of asymmetry and concentration for 
galaxies in Field #1 as function of their visual morphological 
classes. We only show structural values for one field to avoid 
overcrowding in the diagram. Results are similar for Field 
#2. This type of diagram has been extensively used to illus- 
trate and calibrate the morphological segrega tion of gal axies 
based purely on machine class ifications (see Abraha m et alJ 
1996a b; Men anteau etaDll999t) . We use this approach to 
isolate galaxies in the same three equivalent broad visual mor- 
phological categories. The A-C classification recovers <;80% 
of the visual classifications in all three categories, with con- 
taminations of <J15% from other classes. We will use this 
criteria to isolate galaxies based on A-C and compare their 
properties in the next section. 

3. ANALYSIS 

One of the objectives in this paper is to investigate the cos- 
mic star formation history and its dependence on morpho- 
logical type. To accomplish this we need to obtain reliable 
SFRs for a large number of galaxies. Because of the faint 
limits of the sample and the large number of objects involved, 
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FIG. 3. — The star formation rate (SFR) from the spectroscopic sample of 
galaxies derived from the [O II]A3727 EW compared to the SFR estimation 
based on the 1600A UV flux. The inset panel shows the relation between 
the BPZ photometric redshifts estimates and the spectroscopic ones from our 
ground-based observations. 



securing spectroscopic diagnostic for all of them would be 
extremely time consuming, and in most cases unfeasible be- 
cause of their faintness. Instead, we chose to follow an ap- 
proach that takes advantage of the deep and precise ACS mul- 
ticolor photometric information, and we use ancillary ground- 
based spectroscopy and UV imaging for a control sample of 
galaxies to validate our results. 

3.1. The UV SFR at 1600 A 

There are several SFR indicators based on transforma- 
tion between ultraviolet flux and star formation activity 
(seelSteidel et aljH999l: iCoe et al]l200l iHonkins et al]l200lt 
2004, for a discussion and references therein). We opt for 
computing the SFR for each galaxy using the prescription of 
Meureretal. (1999). This is straightforward and computa- 
tionally easy to implement, and derived from a constant SFR 
and a Salteper initial mass function with limits 0.1 - 100 M . 
They define the ultraviolet flux at 1600A, fi600. as the in- 
tegrated light over a square synthetic filter of width 350A 
and central rest-fram e wavelength at 1600A. Based on this, 
Meur er et alJ l[l299) provide a direct transformation between 
a galaxy's SFR and Fi6oo sucn as 

SFR[M yr _1 ] = i -°- 4 ( M i6oo+i8.i5) 

and M1600 is the Fieoo corresponding absolute magnitude in 
AB system. 

Because our ACS filter information only spans between ~ 
3500-9000A, the galaxies' rest-frame UV light only falls into 
our bluest band for z > 1 .5. However, due to the large number 
of objects involved in this study and their precise photometry 
it is possible to derive a general description the galaxies' spec- 
tral energy distribution (SED). We utilize the galaxy's recon- 
structed SED from their BPZ spectral type, 7g, and based on 
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this we estimated the galaxy's UV light at 1600A and compute 
their rest-frame F1600 flux. The 7g-types provided by BPZ 
are a linear combination from a template library of SEDs. We 
use the CW WSB_Benitez2003 template set as described in 
Bemte^ydJ l|2QQ4 ]). These are ba sed o n the templates from 
IColeman. Wu. & Weedmanl (11980ft and iKinnev et all (11996ft 
consisting of El, Sbc, Scd, Im, SB3, and SB2, which represent 
the typical SEDs of elliptical, early/intermediate type spiral, 
late-type spiral, irregular, and two types of starburst galaxies 
respectively. This template set has been modified from ear- 
lier BPZ versions to remove differences between the predicted 
colors and those of real galaxies, whic h results in improved 



BPZ estimates (see Bemte 



for details). We note 



that C oe et al.l d2005) have followed a similar approach to es- 
timate the UV-SFR using also the template fits from BPZ for 
z < 6 galaxies in the HUDF, but at 1400A instead. 

However, because new stars are born in dusty regions, their 
observed UV light can be strongly attenuated and our ob- 
served UV SFRs are therefore subject to obscuration. In order 
to account for this we must correct the SFR estimates for dust 
extinction. We use the empirical relationship between the un- 
obscured far-infrared (FI R) flux and the intrinsic SFR (SFR,) 
prescription by Hop kins et al.l (12001ft . a s FIR light pen etrates 
the dust and provides reliable SFR, . The Hopkins et al. empir- 
ical formulation between the (obscured) observed-UV SFR, 
SFR„, at a given wavelength A and their SFR, is, 



log(SFRf) 



log(SFR o )-0.346x£(A) 

0.797 log(SFR,) + 3.834 



x lo 



(3) 



where SFRs a re in un its of [M©yr '] and k(X) is the empirical 
iCalzetti et all ( 12000ft attenuation law, k(X) =A(X)/E(B-V) 
This improper equation must be solved numerically to ob- 
tain the extinction-corrected SFR, for a given observed UV 
SFR„. We compute the corrected SFR for our full sample for 
A = 1600A using fc(1600) = 9.97. This same method has also 
been used in the HUDF for high redshift galaxies (see Coe 
et al. 2005 and references therein). In the next section we 
use the corrected SFR, to calculate the star formation density 
evolution as a function of morphological type. 

3.2. Emission Line vs. Synthetic SFR 

We compare our SFR estimates obtained using the 1600 A 
flux with independent SFR measurement from traditional 
EWs. For this, we utilize the EW from the [O II]A3727 
emission l i nes to estimate SFR using the prescription from 
iKennicuttl (11998ft using the control-sample of galaxies from 
ourLDDS2 spectroscopic campaign. In Fig|3]we show the in- 
dividual SFR1600 estimates from our synthetic approach and 
the ones obtained from [O II]A3727 EWs. For the majority 
of galaxies there is good agreement between the two mea- 
surements, showing a clear relation between both estimators. 
We find that although there are minor discrepancies for indi- 
vidual galaxies, the overall agreement between SFRigoo and 
[O II]A3727 is quite satisfactory and we estimate out typical 
uncertainties to be ASFR/SFR ~ 0.3. As in this paper we tar- 
get the general properties for the population as a whole rather 
than the precise characterization of each galaxy, we are confi- 
dent that our template-based SFR yield meaningful results. 

3.3. Comparing Synthetic UV fluxes with GALEX 
observations 




f 
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FIG. 4.— The GALEX observed FUV (left panel) and NUV (right panel) 
magnitudes compared to the synthetic recovered values using the recon- 
structed spectral energy distributions from the BPZ template library set. 
GALEX FUV and NUV error bars come from the original errors in the 
GALEX catalogs. Synthetic error bars reflect only original ACS photometric 
errors. 



We further investigate our ability to model the galaxies' 
UV light by comparing our synthetic estimates of the UV 
flux with the observed magnitudes obtained from overlapping 
Galaxy Evolution Explorer (GALEX. iMartin et all2005ft ob- 
servations for the same ACS regions. We retrieved the po- 
sitions, FUV(1530A) and NUV (2310A) magnitudes and as- 
sociated errors in the AB system from the MAST/GALEX 
Web Archive 6 and matched them with the galaxies in our ACS 
sample. From a total of 165 GALEX sources in both fields, 
we successfully identified 71 galaxies in our sample with both 
FUV and NUV magnitudes. We compute the synthetic mag- 
nitudes through the FUV and NUV bandpasses for each of the 
matched galaxies using the reconstructed SEDs. In Figure^ 
we show the comparison of real and synthetic values of FUV 
and NUV. We find that the real and modeled values agree well, 
with a nominal mean scatter A(I/Vgalex — UVsed) of 0.31 
and 0.36 mags for FUV and NUV respectively. Therefore, we 
conclude that our UV flux estimates for ACS galaxies repre- 
sent a good description of their real values and can be used 
as surrogates when measuring the SFR for large samples of 
galaxies. 

4. MORPHOLOGICAL DEMOGRAPHICS BETWEEN 0.3 < Z < 1.2 

4.1. The morphological SFR density 

There is a plethora of studies that have focused on the evo- 
lution of the SFR density, /Q,sfr(V) ? at int e rmediate re dshifts 
llGallego et ail Il995t iMadau et al.l | 1996[ iLillv et alJ 119961: 
iFlores et alJll999HGlazebrpok et alll2004ft and high redshifts 
( Giavalisc o et alJ l2004at iBouwens et alJ 120041 ICoe et all 
2005). These studies have found that SFR rises rapidly from 
z = to z = 1 and peaks somewhere at z > 2 - 3 depending 
on the role of obscuration and our ability to unearth miss- 
ing galaxies at higher redshifts. However, there is little in- 
formation about how the global SF history is distributed into 
galaxy types (e.g. according to morphology). Only a few at- 
tempts have been made to address this issue using relatively 
small, but deep, sample such as the Hubble Deep Fields (see 
Brinchmann et al. 1998; Menanteau et al. 2001). In this paper 
we aim to separate the contributions to the global SFR from 
each of the morphological classes up to z — 1 . Many studies 
have addressed the evolution of psfr(z) in detail and at high 
redshifts, but a full census of the contribution from individual 
morphological classes is yet missing. 

6 http://galex.stsci.edu/GRl 
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FIG . 5 . — The co-moving star formation density of galaxies as a function of redshift for the three morphological groups. Red rounded squares represent E/SOs, 
green squares, Spirals and blue circles Irr/Pec selected galaxies. Bars represent 1<t deviations. Open triangles represent the total star formation density from all 
classes. The left panel shows the results from morphologically selecting galaxies according to their visual classes, while the right panel shows the results from 
the A-C automatic selection. 



We calculate the SFR density for each morphological class 

at given redshift interval z\ < z < Z2 as, 

, , _ SFR, 
Psfr(z)= £ — - 1 (4) 

j,Zl<Zj<Z2 VrmX 'j 

where we sum the star formation contribution of each galaxy, 
SFR 7 , for zi < Zj < Z2 and V max is the co-moving volume 
at Zmax, the highest redshift at which a given galaxy, j, is 
still brighter than our sample limiting magnitude, zsso.ab < 25 
iSchmidtl[l9 68). We proceed in this fashion for all galax- 
ies, and calculate psfr(z) f° r the E/SO, Spirals and Irregu- 
lar/Peculiar samples. We explore the morphological selection 
calculating psfr(z) using both the visual morphological clas- 
sifications and the machine-based selection. In Fig.[5]we show 
the results of ,osfr(z) for the visual and automatic classifica- 
tions. We also show the total psfr(z) as well as the contri- 
bution of the three subclasses. We note that E/SOs are, as 
expected, the lowest contributors to the total Psfr(z) regard- 
less of the classification used. We see from the figure that 
there are only small differences in psfr(z) from the two ap- 
proaches, mostly in the selection of E/SO and Irr/Pec. More- 
over, the star formation density for early-types shows a mod- 
est increase with redshift which seems to flatten out by z — 1 . 
This is not the case of either of the other classes or the global 
Psfr(z)- On the other hand the Irr/Pec galaxies show a con- 
stant and in step increase with redshift virtually matching the 
contribution from Spirals, the most active group, by z — 1 . 

4.2. Differential Evolution 

We investigate the morphological evolution of galaxies us- 
ing the observed color magnitude diagrams at several red- 
shift intervals. Our motivation arises from the recognition 
that early-types in clusters exhibit a tight correlation be- 
tween their colors and magnitudes. This correspondence, 
known as the color-magnitude relation (CMR) or "red se- 
quence", has been widely used to constrai n the ages of galax- 
ies, particularly for cl uster ellipticals (see lBower et alJ[l992t 
iBlakeslee et alJl2~003bl and references), and in re cent years 
also for field galaxies at different cosmic times JBell et alJ 
l2004lKodama. Bower. & Bell 1999ft . Galaxies at similar evo- 
lutionary stages should have comparable colors, while depar- 
tures from this can be employed to establish variations in their 



star formation history and be interpreted as differential evolu- 
tion within their stellar population, regardless of its origin (i.e 
secular evolution or recent mergers). 

We construct the color magnitude (CM) diagrams for (775 - 
Z850 and V(,o6 ~ hi5 colors as a function of z&sq at three evenly 
spaced redshift intervals, 0.3 < z < 0.6, 0.6 < z < 0.9 and 
0.9 < z < 1.2. We chose to avoid the z < 0.3 objects because 
they are highly incomplete and undersampled. Each interval's 
diagram is a snapshot of the galaxies' current evolutionary 
phase, where dispersions in their color can be interpreted as 
variations (or lack of) in their age and star formation history. 
Despite the relative large width of our redshift bins, we can 
use them to provide a general census of the galaxies' relative 
formation stages as a function of their morphology. 

In Fig [6] we plot the diagrams keyed to their visual mor- 
phological type, for V(,os~hi5 m the left panel and (775-2850 
on the right one. We use the same symbols as in Fig. |2] 
When focusing on the red envelope of galaxies, we see that 
this is dominated by morphologically selected E/SO galax- 
ies in all redshift bins. Although there is considerable scat- 
ter, a red sequence of passively evolved E/SO is quite distin- 
guishable. We also show the z = 0.023 Coma Cluster CMR 
( Bowe r et al J[l9 92 ) transformed to the intervals' central red- 
shift and transformed from U — V to (775-2850(2) and V606 - 
his(z) accordingly. We use empirical and synthetic tem- 
plates dBruzual & Charlotl 120031 IColeman. Wu. & Weedmanl 
1980) to transform between A(U -V) to A(i775 -Z85o) and 
A(V606-z'775) respectively. The redshifted Coma CMR rep- 
resents a bona-fide benchmark for old and coeval early-type 
systems. For the highest redshift (775 -Z850 interval we plot 
the empiri cal CMR for the z = 1 .23 7 CL1252-2927 cluster as 
derived by Bla keslee et alJ J2003 b). In each of the panels we 
also show with a solid st ar the reds hifted values of Mg from 
the COMBO 17 survey JWolf et al.l2 003) to z 8 50 at the central 
redshift of each interval. This insures that even in our highest 
redshift bin we are probing galaxies at ~ 2 magnitudes deeper 
than M* within the flux limit of our sample. 

Previous CMR studies of HST morphologically 
selected field galaxie s show similar results. The 
Kodama, Bower, & Bell (1999) analysis of ellipticals in 
the Hubble Deep Field North at (z) ~ 0.9 was the first to 
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FIG. 6. — The V606 _ '775 and '775 _ Z850 observed color magnitude diagrams as a function of Z850.AB magnitude for three evenly spaced redshift intervals between 
Z = 0.3- 1.2. Symbols are the same as in Fig.|2|and are keyed to their visual morphological classes. E/SO systems are represented by ellipses; Spiral galaxies 
by spirals and Irr/Pec as asterisks. The dashed lines represent the color magnitude relation for E/SOs for t he Coma Cl uster redshifted to mean interval redshift of 
each panel. The solid line in the right panel is the empirical relation for CL1252-2927 at z = 1.237 from Blakeslee et al. 1 2003b). For comparison we show as 
solid stars the values of M* B from the COMB017 survey I Wolf et al. 2003) for early types redshifted to z&so at the mean redshift of each interval. 



report the presence of a red sequence of field ellipticals at 
that redshift, comprising around ~ 1/2 of the population 
of early types. However, a significan t fraction have bluer 
rest-frame colors, which Menanteau et al. (2001) later re- 
ported to be also coincidental with havin g large internal 
color dispersions. In a later larger sample, Bell et al. (2004) 
using GEMS observations confirmed the presence of a red 
sequence dominated by early types (70% -80%) at z ~ 0.7. 
These studies establish a pattern in morphologically selected 
samples of spheroids where systems are predominantly old 
and coeval, but bluer and presumably younger ellipticals 
coexist at i ntermediate redshifts (see Menanteau et al.j|1999l 
[Kodama et alJ 119991 lAbraham et al.1 Il999t iM enantea u et all 
l200lU2004HBenson et all2002HCross et all2004l) . 

In this paper we report a similar behavior. We note that 



while E/SOs tend to dominate the red envelope, there is an 
important fraction of systems with significantly bluer colors 
than those of old ellipticals. Moreover, redder E/SOs tend 
to be brighter (massive) than bluer ones, which are system- 
atically fainter at all redshifts. This result favors a view in 
which most massive ellipticals might have formed at higher 
redshift and where at least an important fraction of those were 
already in place at z — 1. However, it is unclear whether 
bluer early-type galaxies will become the present day mas- 
sive systems we see in place. We note that it is unlikely 
that these bluer systems could be confused with the bright- 
est knots of faint irregular galaxies, as it has been established 
that ACS can efficiently distinguish between such contrast- 
ing galaxy types at these li mits, and especial ly in these high 
signal-to-noise images (see lCross et al.l2004l:IConsehce et alJ 
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l2004t lMenante a~et al.l20 04). Spiral galaxies display a broad 
range of colors, and nearly uniform distribution of magnitudes 
overall redshift ranges with observed magnitudes similar to 
those of the brightest E/SOs in the red envelope. Irregular and 
Peculiar galaxies, on the other hand, tend to dominate the low 
luminosity loci with predominantly bluer colors, but compa- 
rable to the blue E/SO systems. While these two groups — 
blue E/SO and Irregular/Peculiars — have dramatically differ- 
ent structural properties, their similarities in colors and lumi- 
nosities make it tantalizing to relate the evolutionary histories 
of these two very different types. 

5. SUMMARY 

In this paper we have investigated the properties of a 
large morphologically selected sample of ACS galaxies. The 
deep multicolor information and ancillary ground-based spec- 
troscopy has enabled us to constrain the evolution of the SFR 
density and the individual contribution to the global psfr(z)- 
We report that spiral galaxies are the main providers to the 
global SFR density at all redshifts, while the E/SOs' humble 
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